14th  International  Conference  on  Information  Fusion 
Chicago,  Illinois,  USA,  July  5-8,  2011 


Optimal  Placement  of  Heterogeneous  Sensors  in 

Target  Tracking 


Chun  Yang 

Sigtem  Technology,  Inc. 
San  Mateo,  CA  94402 
chunynag@sigtem.com 


Lance  Kaplan 

ARL/SEDD 
Adelphi,  MD  20783 
lance.m.kaplan@us.anny.mil 


Erik  Blasch 

AFRL/RYAA 
WPAFB,  OH  45433 
erik.blasch@wpafb.af.mil 


Michael  Bakich 

AFRL/RYAA 
WPAFB,  OH  45433 
michael.bakich@wpafb.af.mil 


Abstract  -  Target  tracking  performance  is  determined  by 
the  fidelity  of  target  mobility  model  (F,  Q),  tracking 
sensor  measurement  quality  (R),  and  sensor-to-target 
geometry  (H).  A  tracking  sensor  manager  has  choices  in 
sensor  selection/placement  (H),  waveform  design  (R),  and 
filter  tuning  (F  and  Q),  thus  affecting  the  tracking 
performance  in  many  ways.  This  paper  concerns  with  the 
geometry  aspect  of  sensor  placement  so  as  to  optimize  the 
tracking  performance.  Recently,  a  considerable  amount 
of  work  has  been  published  on  optimal  conditions  for 
instantaneous  placement  of  homogeneous  sensors  (same 
type  and  same  measurement  quality)  in  which  the  targets 
are  either  assumed  perfectly  known  or  the  target  location 
uncertainty  is  averaged  out  via  the  expected  value  of  the 
determinant  of  the  Fisher  information  matrix.  In  this 
paper,  we  derive  conditions  for  optimal  placement  of 
heterogeneous  sensors  based  on  maximization  of  the 
updated  Fisher  information  matrix  from  an  arbitrary 
prior  characterizing  the  uncertainty  about  the  initial 
target  location.  The  heterogeneous  sensors  can  be  of  the 
same  or  different  types  (ranging  sensors,  bearing-only 
sensors,  or  both).  The  sensors  can  also  make,  over  several 
time  steps,  multiple  independent  measurements  of 
different  qualities. 

Keywords:  Performance  Evaluation,  Fisher  Information 
Matrix,  Sensor  Placement,  Heterogeneous  Sensors, 
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1.  Introduction 

Target  tracking  performance  is  determined  by  the  fidelity 
of  target  mobility  model  (F,  Q),  tracking  sensor 

measurement  quality  (R),  and  sensor-to-target  geometry 
(H).  The  effectiveness  of  a  sensor  update  is  implicitly 
dependent  of  the  range  to  target  via  signal  to  noise  ratio 
(SNR)  for  a  ranging  sensor  while  explicit  for  a  bearing- 
only  sensor.  A  tracking  sensor  manager  has  choices  in 
sensor  placement  (H),  waveform  design  (R),  and  filter 
tuning  (F  and  Q)  to  control  tracking  performance.  This 
paper  is  concerned  with  the  geometrical  aspect  of  sensor 
placement  so  as  to  optimize  the  tracking  performance. 

Placement  of  m  sensors  around  a  target  has  drawn  a 
considerable  amount  of  attention  recently.  Interesting 
results  are  summarized  in  a  recent  paper  [1]  where  the 


determinant  of  the  Fisher  information  matrix  (FIM)  is 
maximized  so  as  to  obtain  necessary  and  sufficient 
conditions  for  optimal  placement  of  ranging  sensors, 
bearing-only  sensors,  and  time  of  arrival  (TOA)  and  time 
difference  of  arrival  (TDOA)  sensors,  respectively. 
Similarly,  the  use  of  the  trace  of  the  Cramer-Rao  lower 
bound  (CRLB),  which  is  the  inverse  of  the  FIM,  is 
considered  in  [5,  6,  7]  for  TDOA  measurements. 

As  pointed  out  in  [2],  a  majority  of  work  in  optimal  sensor 
placement  including  the  references  cited  above  make  a 
fundamental  assumption  that  the  target  location  is  known 
perfectly,  which  is  hardly  practical  but  establishes  some 
theoretical  bounds.  When  the  uncertainty  in  target  location 
is  characterized  by  a  probability  distribution  such  as 
truncated  radially  symmetric  distributions,  the  expected 
value  of  the  FIM  determinant  is  used  in  [2]  to  obtain  useful 
results  with  a  combination  of  analysis  and  simulation. 

Furthermore,  the  above  references  assume  homogeneous 
sensors.  That  is,  they  are  of  the  same  type  and  of  the  same 
quality  for  their  respective  necessary  and  sufficient 
conditions  of  optimal  placement  to  hold. 

In  this  paper,  we  consider  the  problem  of  optimal 
placement  of  heterogeneous  sensors  in  tracking  of 
uncertain  targets.  By  heterogeneous,  we  mean  the  use  of 
sensors  that  are  different  in  sensor  type  and  measurement 
quality.  In  addition,  to  account  for  target  uncertainty,  we 
use  the  updated  FIM  to  derive  optimal  placement 
conditions. 

Different  from  previous  work,  which  is  mostly  formulated 
in  a  static  snapshot  fashion,  we  consider  a  tracking 
application  in  which  the  time  factor  is  involved  and  the 
sensor  placement  actually  occurs  in  time  and  in  space.  As 
an  example,  a  number  of  UAVs  receive  a  handover 
message  from  a  surveillance  aircraft  in  a  layered  sensing 
scenario  and  they  coordinate  to  track  a  designated  target  in 
a  cooperative  manner.  Since  a  UAV  is  flying  at  a  constant 
altitude  and  solves  for  a  ground  target  location,  a  2-D 
model  is  reasonable  and  will  be  considered  as  in  most 
sensor  placement  literature. 

Assume  that  the  tracking  sensors  under  consideration  will 
make  a  unit  movement  per  unit  time  and  each  movement 
has  an  associated  cost.  Due  to  higher  risk  for  covertness 
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and  survivability,  it  is  reasonable  to  assume  that  a  radial 
movement  to  get  closer  to  a  target  has  a  higher  cost  while 
tangential  movement  carries  a  constant  cost.  A  UAV  thus 
faces  four  choices: 

•  Do  not  move  and  take  independent  measurements.  If  the 
current  geometry  is  favorable,  the  accumulation  of 
independent  measurements  has  the  effect  of  reducing 
measurement  errors.  A  steady  state  is  soon  reached  and 
no  further  reduction  can  be  expected  for  this  geometry. 

•  Move  radially  toward  the  target.  This  runs  a  higher  cost 
for  movement  but  has  the  benefits  of  increasing  SNR  and 
improving  the  measurement  quality,  particularly  for 
bearing-only  sensors. 

•  Move  tangentially  around  the  target.  This  allows 
maximum  change  in  viewing  geometry  if  moving  in  the 
direction  of  the  principal  axis  of  the  error  ellipse. 

•  A  combination  of  all  the  above. 

The  above  approach  decides  a  move  from  one  time  step  to 
next  and  is  therefore  a  greedy  one.  The  objective  is  to  seek 
the  best  outcome  in  terms  of  best  accuracy  in  minimum 
time  and  cost  with  least  risk. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2, 
the  sensor  models,  optimality  criteria,  and  information 
updating  equations  are  introduced.  Section  3  presents  three 
strategies  for  placing  two  homogenous  and  heterogeneous 
sensors,  either  co-located  or  distributed.  Simulation  results 
are  analyzed  in  Section  4.  Finally,  concluding  remarks  are 
provided  in  Section  5. 

2.  Information  Updating 

Consider  a  target  with  an  unknown  state  ( e.g .,  position  and 
velocity)  x  and  the  i-th  sensor  with  a  known  state  x,.  The  i- 
th  sensor’s  measurement  is  given  by 

A=/''(x>x/)+v,  (!) 

where  fl[  (■,•)  is  typically  a  nonlinear  equation,  the 

superscript  /,-  e  {r,  6,  <p,  ...}  is  a  label  indicating  a 
particular  type  of  nonlinearity  for  the  ;-th  sensor:  range, 
azimuth,  and  elevation,  and  v,-  is  the  sensor  measurement 
error,  assumed  to  be  a  zero-mean  Gaussian  with  variance  ( 
a')2,  denoted  by  27sT{0,  (a'1)2}.  In  this  paper,  the 
dimension  of  the  state  vector  x  is  n. 

Given  m  such  sensors,  the  relationship  between  the 
measurements  and  target  state  x  can  be  represented  in 


matrix-vector  form  as 

z  =  f(x)  +  v  (2) 

where 

z7=[z,  z2  •••  zm\  (3a) 

vr=[v,  V2  ...  vj  (3b) 

fO)  =  |/',(x,x1)  /'"(x,x„,)J  (3c) 


and  v  ~  .Afi).  R)  where  the  superscript  T  stands  for  vector 
or  matrix  transpose. 

For  an  unbiased  estimate  x  of  x,  the  Cramer-Rao  lower 
bound  (CRLB)  states  that: 

£{(x  -  x)(x -  x)r}  >  CRLB  =  M~l (x)  =  T(x)  (4a) 

where  'M(x)  is  called  the  Fisher  information  matrix  (FIM) 
and  its  ijth  element  is  defined  as: 

[2M(x)]iy  =  (x)f  ^-ln(f  (x))}  (4b) 

dbc,.  dXj 

It  is  clear  from  (4a)  that  if  JM(x)  is  singular,  no  unbiased 
estimator  exists  for  x  with  a  finite  variance.  When  the 
measurement  noise  is  Gaussian,  the  FIM  is  given  by: 

M(x)  =  VrJ(x)R'VJ(x)  =  UTR-'H  (5) 

where  Vxf  is  the  Jacobian  of  the  measurement  vector  f  with 
respect  to  x,  denoted  by  H. 

In  [1],  the  determinant  of  the  FIM  is  maximized  so  as  to 
obtain  necessary  and  sufficient  conditions  for  optimal 
placement  of  ranging  sensors,  bearing-only  sensors,  and 
time  of  arrival  (TOA)  and  time  difference  of  arrival 
(TDOA)  sensors,  respectively.  Since  f  is  also  a  function  of 
{x,-,  ;  =  1,...,  m},  the  optimization  can  be  formulated  as: 

{xiv,xI,}  =  arg  max  det(jM(x,x1,...,xm)) 

t*i . x.) 

=  arg  max  d ct (l  I '  R1 1 1  )  (6) 

t*i. 

Note  that  the  FIM  is  evaluated  at  x,  which  implicitly 
assumes  that  the  target  location  is  known  perfectly.  In 
practice,  however,  only  a  coarse  estimate  of  the  target  state 
is  available.  Assume  that  the  prior  on  the  target  state  is 
characterized  by  a  distributed  denoted  by  /j(x).  The 
expected  value  of  the  determinant  of  the  FIM  is  used  in  [2] 
for  optimal  sensor  placement: 

{x*,...,x* }  =  arg  max  }  J  det(jM(x,  Xj,...,  xm))p(x)dx  (7) 

In  this  paper,  we  consider  a  Gaussian-distributed  initial 
target  state  as  x  ~  dA(x0,  P0)  and  perform  optimal  sensor 
placement  based  on  the  updated  FIM  evaluated  around  the 
predicted  state  x0. 

Given  the  initial  target  state  distribution  as  x  ~  p(x)  = 
AVfxo.  Po),  the  maximum  a  posterior  (MAP)  estimator  is 
actually  a  nonlinear  least  squares  estimator  given  by 

x  =  argmin[(z-f(x))rR_l(z-f(x))+(x-x0)rP0_1(x-x0)J  (8) 

X 

The  covariance  of  the  state  estimate  can  be  approximated 
by  linearizing  f(x)  about  the  mean  of  the  prior.  Specifically, 
one  assumes  that 

f(x)  =f(x0)  +  H(x  —  x0)  (9a) 

Hr  =  [h,  hj  (9b) 
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d/'‘ 

\% 
_ 1 

lx0 

5xj 

dx2 

5x„ 

(9c) 


where  V./ stands  for  the  gradient  of  the  scalar  function/ 
with  respect  to  x. 

By  inserting  (9)  into  (8),  it  is  easy  to  show  that  the  linear 
least  squares  estimate  of  the  target  state  is 

x  =  x0  +  (HrR  ]H  +  P/r'H7 R-‘(z  -  f(x0))  (10) 

and  the  covariance  of  the  state  estimate  for  the  posterior 
distribution  is 

P  =  £{(x-x)(x-x)r}  =  (HrR"1H  +  P0“1)“l  (11) 

The  covariance  given  by  (11)  is  also  the  CRLB  derived 
from  the  joint  distribution  of  the  target  state  and 
measurements.  From  (4a),  M  =  P"1  and  the  information 
updating  equation  is  given  by 

M  =  M0+HrRH.  (12) 

Usually,  the  sensors  are  spatially  disparate  so  that  the 
measurement  errors  v,  are  statistically  independent  with 
variance  cr2 .  Then,  the  information  update  equation  can  be 

rewritten  as 


M  =  M0  + 


=M0+jr4Ah  j 


where 


h 

_ i_ 

llh,.  I 


and  ^2 


(13a) 


(13b) 


Inspection  of  (13)  reveals  that  each  sensor  provides  a  rank- 
one  positive  definite  addition  to  the  information  matrix, 
which  increases  the  “size”  of  M.  The  unit-norm  £ 

represents  the  angular  orientation  of  the  single  sensor 
contribution  in  information  space,  and  the  reciprocal  of  fff 

represents  the  power  (or  quality)  of  the  sensor 
measurement. 


Table  1  provides  the  value  of  £  and  ff2  for  2-D  range  and 

bearing  measurements  in  Cartesian,  polar  (r,  9),  and 
spherical  (r,  0,  </>)  coordinates,  respectively.  In  the  2-D 
case,  the  angle  6  of  the  LOS  from  the  sensor  to  target  is 
measured  relative  to  the  x-axis.  Table  1  reveals  that  the 
orientation  £  of  range  sensors  align  with  the  LOS  vector 

between  the  sensor  and  the  target.  On  the  other  hand,  the 
orientations  of  bearing  sensors  are  orthogonal  to  this  LOS 
vector.  The  quality  of  the  range  sensors  only  depends  on 
the  actual  measurement  variance.  However,  for  a  given 
angular  measurement  error,  the  quality  of  sensor  degrades 
as  the  distance  to  the  target  increases  because  the 
corresponding  state  value  is  a  cross-range  quantity. 
Furthermore,  it  is  reasonable  to  assume  that  the 
measurement  error  also  increases  as  the  target/sensor 
distance  expands  due  to  a  decrease  in  the  SNR.  Without 


accounting  for  environmental  effects,  the  measurement 
error  for  a  point  target  is  ideally  invariant  to  the  angular 
coordinates.  As  a  first  order  approximation,  it  is  reasonable 
to  assume  that  £  and  a1  only  depend  on  the  angular  and 

range  coordinates  of  the  sensor,  respectively,  in  a  polar 
coordinate  system. 

Table  1 .  Sensor  Measurement  and  Linearized  Equations  in  2D 


(a)  2D  Active  Sensors  with  Ranging  Measurements 

/Hx,xi)  =  1/(x-xi)3+o-^)s  a?  =  of 


k,= 


yg-yj 


=  [cos(q)  sin(3)f 

n  =  n/Cxo-TTkSvT)7 

()>)  2D  Passive  Sensors  with  Bearing-Only  Measurements 

=  [-  sin(  Q)  cos(6j)]r 


3.  Optimal  Placement  in  Time  and  Space 

For  2D  scenarios,  a  minimum  of  two  sensors  are  sufficient 
to  obtain  an  optimal  solution.  The  two  sensors  can  be  of  the 
same  type,  observing  a  target  from  two  different  directions 
at  the  same  time.  Or  the  two  sensors  can  be  of  different 
types,  observing  the  target  from  the  same  direction  (co¬ 
located)  at  the  same  time.  This  is  an  instantaneous 
placement  with  single-look  optimality.  For  stationary  or 
slow-moving  targets,  sensors  can  stay  at  the  same  place  or 
move  to  different  locations  while  taking  multiple 
measurements.  This  is  a  temporal  placement  with  multiple- 
look  optimality.  Conditions  for  optimal  placement  are 
derived  below. 

3.1  Optimality  Criteria 

Assume  that  the  number  of  sensor  measurements  is 
constrained  to  be  m.  From  (13a),  it  is  easy  to  see  that  the 
trace  of  the  updated  information  matrix  for  homogeneous 
measurements  is 


T  =  trace(M0)  +  -^- 
a 

This  trace  of  the  information  matrix  accumulates  the  sensor 
quality  of  all  sensor  updates  involved  without  regard  to  the 
effects  of  angular  orientation  of  the  sensors.  The  following 
theorem  places  a  lower  bound  on  the  achievable  updated 
covariance  matrix  P. 

Theorem  1.  Consider  an  «x«  covariance  matrix  P.  Let  M  = 
P'1  be  the  information  matrix.  If  traced M)  =  T  is  constant, 
then  trace( P)  >  n2IT,  det( P)  >  ri'IT'.  The  equality  occurs  if 
M  =  (77/7)1  and  P  =  {n/T)\. 

In  short.  Theorem  1  states  that  if  the  updated  covariance 
matrix  has  equal  eigenvalues,  no  other  configuration  of 
sensors  could  provide  a  tighter  covariance  error  when  the 
trace  of  the  information  matrix  is  fixed.  In  this  case,  the 
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error  ellipse  is  actually  circular  for  the  2D  cases.  Theorem 
1  is  a  generalization  of  a  theorem  in  [3]  and  its  proof  is 
given  in  [8], 


circular.  The  condition  indicates  that  two  placements  such 
that  0\  -  &i  =  90°  provide  an  optimal  solution,  which  is 
consistent  with  our  intuition. 


In  general,  each  measurement  can  lower  one  eigenvalue  of 
Po-  For  an  /7-dimensional  state  vector,  it  might  be  possible 
for  7?-l  measurements  to  lower  n- 1  eigenvalues  to  the  value 
of  the  lowest  X\.  Flowever,  depending  on  the  spread  of  the 
eigenvalues,  more  measurements  might  be  necessary  to 
achieve  a  desired  spherical  error.  Even  if  sensor  placement 
cannot  achieve  a  spherical  state  error,  Theorem  1  provides 
intuition  that  it  is  best  to  try  to  achieve  as  close  to  uniform 
eigenvalues  in  the  updated  covariance  as  possible. 

3.2  Two  Homogeneous  Sensors 


Bringing  the  first  condition  (16b)  into  (15b)  gives: 

(17a) 
(17b) 

That  is,  the  optimal  placement  is  for  the  two  sensors  at  0X  = 
-02  (symmetric  about  a  principal  axis  of  the  error  ellipse). 
The  condition  for  such  an  optimal  placement  is: 


sin'  9,  =  - 


1  f  o2(a-b ) 
2k 


+  1 


#,=-#,  =  sin  I— 


->  1  a2(a-b) 


2k 


+  1 


In  this  section,  we  consider  two  homogeneous  sensors  that 
are  of  the  same  type  and  same  quality.  We  will  use  the 
unified  sensor  model  (13b)  that  accounts  for  both  ranging 
and  bearing-only  measurements  as  listed  in  Table  1(a)  and 
(b).  Assume  that  the  prior  on  the  target  state  is  x  ~  AVfx(l, 
Po).  Around  x0,  we  can  linearize  the  measurement 
equations  for  given  xj  and  x2,  leading  to  observation 
matrices  h,  and  hb  respectively.  By  a  proper  coordinate 
rotation,  we  can  diagonalize  P0  such  that  we  have  a 
diagonal  information  matrix  prior  M0  =  Po  1  =  diag([a,  6]). 

Given  the  prior  M0  =  diag([a,  ft]),  we  want  to  place  two 
sensors  at  h  =  [cos#i,  sin#i]T  and  h,  =  [cos#>,  sin#>]T 
with  an  equal  effective  measurement  error  variance  a2  and 
to  take  a  number  of  k  independent  measurements  in  the 
placement  so  as  to  accumulate  the  needed  gain  (accuracy). 
What  are  k,  0X,  and  0,  that  result  in  the  maximum  updated 
information  matrix? 
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It  is  easy  to  verify  that  when  a  =  b,  the  angular  sector 
between  the  sensors  is: 


0,  -02  =2sin-‘^  =  90°  09) 

which  is  consistent  with  the  second  condition  (16b). 

When  |  a2  (a  -  b)\  >  2,  there  is  no  single  look  update  (i.e.,  k 
=  1)  to  achieve  an  instantaneous  optimality.  However,  a 
number  of  independent  updates  can  be  used  to  accumulate 
the  required  gain  to  ensure  a  solution.  The  number  is  given 
by: 


According  to  (13a),  the  updated  information  matrix  can  be 
written  as: 


M  =  M0  +  4-(h,£f  +  h,h2r)>£>1 

a 

k  k 

a  +  —  (cos2  0X  +  cos2  02 )  -737- (cos  6X  sin  6X  +  cos  02  sin  02 ) 

_  G  <7 

k  k  2  2 

—r  (cos  6X  sin  6X  +  cos  02  sin  02 )  b  +  —  (sin  0X  +  sin  02 ) 
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(14b) 


To  maximize  the  updated  information  matrix  is  equivalent 
to  making  (14b)  a  scaled  identity  matrix.  To  do  so,  we  need 
to  find  k,  6\,  and  0  such  that: 

cos  0  sin  9X  +  cos  02  sin  0  =  0  (15a) 

a  +  -^-(cos2  9,  +  cos"  #,)  =  #  +  —-(: sin2  9,  +sin2  9 1)  (15b) 

cr  o 


(20) 

where  sup{-}  is  the  smallest  of  all  integers  that  satisfy  the 
condition  and  [•]  stands  for  the  least  integer. 

When  the  condition  for  a  solution  (18b)  is  met,  the  angular 
placement  9\  varies  as  a  function  of  k.  In  the  limit  k  goes  to 
infinity,  the  placement  solution  is  the  same  as  (19)  with  9X 
=  -0  =  45°.  That  is,  the  two  sensors  maintain  an  angular 
separation  of  90°  and  symmetric  about  a  principal  axis  of 
the  error  ellipse. 

The  above  analysis  leads  to  an  optimal  strategy  to  the 
placement  of  two  cooperative  sensors.  From  (15b),  we 
have: 


k *  =sup<^  k  > 
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From  the  condition  (15a),  we  have: 

sin?^  +  sin 2#,  =  sin(#,  +  #2)cos(#!  -  02)  =  0  (16a) 

f#,+#2=0  (16b) 

}  #,  -02  =  90” 

The  condition  in  the  second  row  of  (16b)  is  not  valid  in 
view  of  (15b)  unless  a  =  b.  When  a  =  b,  the  prior  is 


cos2  #,  -  sin2  #[  +  cos2  #,  -  sin2  #,  =  — — a^a  (21a) 

k 

cos2#,  +cos2#,  =  — — (21b) 
k 

Adding  the  sums  of  (16a)  squared  and  (21b)  squared  gives: 
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cos  2#,  cos  2#,  +  sin  2#,  sin 2(9,  _  ^  a)  a —  j 
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From  the  first  condition  of  (16b),  0\  +  #>  =  0,  we  finally 
have: 
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The  condition  that 
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leads  to  the  same  condition  as  (20).  It  can  be  shown  that 
(2 If)  and  (17b)  are  equivalent. 

Note  that  the  optimal  placement  condition  derived  in  this 
section  is  for  a  generalized  sensor  with  h  and  a2 ,  which 
can  represent  either  a  ranging  error  or  a  bearing-only 
sensor.  In  the  latter  case,  if  a  sensor’s  LOS  vector  to  target 
has  an  angle  <p  with  respect  to  the  x-axis,  then  the 
observation  vector’s  angle  is  given  by  0=  <p  ±  n  12.  The 
equivalent  measurement  error  is  given  by  a2  =  rcr2  as  listed 
in  Table  1. 


The  determinant  of  M  is  given  by: 
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The  goal  is  to  maximize  det(M)  (equivalent  to  minimize 
det(P)).  Since  the  second  term  on  the  right  hand  side  of 
(24b)  is  always  negative,  it  can  be  nullified  when  0=  0°  or 
90°. 


When  0=  0°  or  90°,  (24b)  becomes: 


det(M)  = 


k  )  k 
a  +  —  b  + 
a,  r  cr. 
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a  +  -^Ab  +  ??\  9  =  9  0 
07 


The  difference  between  the  two  values  of  (25)  can  be 
evaluated,  for  a  >  b,  as: 


3.3  Two  Co-Located  Heterogeneous  Sensors 

In  this  section,  we  consider  two  sensors  co-located  on  a 
same  sensor  platform.  The  two  sensors  are  of  different 
types  with  different  measurement  qualities,  one  providing 
ranging  and  the  other  bearing-only  measurements.  Given  a 
prior  about  the  unknown  target  x  ~  2A(\ 0,  P0),  we  want  to 
determine  where  to  place  the  sensor  platform  x  [  relative  to 
x0  and  to  take  how  many  independent  measurements  so  as 
to  obtain  maximum  updated  information  matrix. 

Again,  the  coordinate  axes  can  be  rotated  to  align  with  the 
principal  axes  of  Po,  leading  to  the  prior  in  the  diagonal 
form  denoted  by  M0  =  Po"1  =  diag([a,  &]).  Around  x0,  we 
can  linearize  the  measurement  equations  for  the  two 
sensors,  leading  to  observation  matrices  h  =  [cos#,  sin#j' 
and  h2  =  [sin#,  cos#|r  with  measurement  errors  0\  and  a?, 
respectively.  What  are  k  and  #  that  result  in  the  maximum 
updated  information  matrix? 

To  start,  the  updated  information  matrix  is: 


M  =  M0  -1-  — V  h , h ; 
07 


=rrh2h2 
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where  r  =  ||X]  -  x0||2  is  the  predicted  range  from  the  sensor 
to  target.  Eq.  (22)  can  be  further  written  as: 
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>  0,  r2<j2  <  <T,2 

=>  #  =  0 

<  0,  r2<j\  >  of 

=>  #  =  90 

(26b) 


The  above  analysis  indicates  that  the  number  of  updates  k 
does  not  affect  the  optimal  placement  and  the  optimal 
sensor  pointing  is  either  along  the  first  principal  axis  or 
perpendicular  to  it.  When  a  >  b,  that  is,  the  largest 
eigenvalue  along  the  y-axis,  the  sensor  should  be  placed  at 
#=  0°  if  rrr22  <  (T-r.  that  is,  the  cross-ranging  error  (along 
the  j-axis)  is  smaller  than  the  ranging  error  (along  the  x- 
axis).  This  strategy  advocates  updating  the  largest  error 
direction  with  the  best  measurement.  In  essence,  sensors 
are  added  such  that  the  posterior  M  remains  diagonal  and 
the  eigenvalue  spread  is  as  small  as  possible. 

3.4  Two  Heterogeneous  Sensors 

In  this  section,  we  consider  a  more  general  case  with  two 
heterogeneous  sensors.  We  start  with  diagonal  information 
matrix  prior  M0  =  diag([a,  #]),  which  is  obtained  by 
rotating  the  coordinate  axes  along  with  the  principal  axes  of 
an  arbitrary  prior. 

Two  sensors  are  placed  at  £j  =  [cos#,,  sin#,]T  and  = 

[cos#,,  sin#>]T  with  their  effective  measurement  errors 
being  ff2  and  ff2 ,  respectively.  Assume  a  number  of  k 
independent  updates  are  taken  in  the  placement  to 
accumulate  the  needed  gain  (accuracy). 
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The  resulting  updated  information  matrix  can  be  written  as: 


M  =  M0+£(4rhX  +4rh,h[)’^  1 

<7,  (7. 
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To  make  (27b)  a  scaled  identity  matrix,  we  need  find  k,  ft, 
and  ft  such  that: 


£cosftsinft  icosftsinft 
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From  the  condition  (28a),  we  have: 
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From  (28b),  we  have: 
cos2ft  cos2 9,  _b-a 


OT 


ft 


cos' 2ft  2cos2ftcos2ft  cos' 2ft 

^4  1  ~2  ~2  1  ^4 


)  —  a 
k 


(30a) 

(30b) 


The  sum  of  (29b)  and  (30b),  after  easy  manipulation,  leads 
to: 


single  look  update  (i.e.,  k  =  1)  to  achieve  an  instantaneous 
optimality. 

But  a  number  of  independent  updates  can  be  used  to 
accumulate  the  required  gain  to  ensure  a  solution.  From 
(31a),  the  number  is  given  by: 
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The  upper  bound  in  (34)  only  limits  k  to  achieve  a  scaled 
identity  information  matrix,  not  on  the  ability  to  further 
increase  the  information  matrix,  which  has  no  limit  on  k. 

Within  the  bounds  in  (34),  for  each  k,  there  is  a  pair  of 
optimal  angular  placements  {ft,  ft}  as  given  by  (32d)  and 
(32e). 

Since  the  optimal  placement  condition  is  derived  for 
generalized  sensors,  h  and  cf2  can  represent  either  a 
ranging  error  or  a  bearing-only  sensor.  Again  in  the  latter 
case,  if  a  sensor’s  LOS  vector  to  target  has  an  angle  ^>with 
respect  to  the  x-axis,  then  the  observation  vector’s  angle  is 
given  by  9=  (p±  n!2.  The  equivalent  measurement  error  is 
given  by  ff2  =  ro1  as  listed  in  Table  1. 

3.5  Generalization  of  Existence  Conditions 

In  [1],  a  condition  on  ranges  of  bearing-only  sensors  is 
formulated  (see  (43)  in  Theorem  9,  [1]),  which  determines 
if  it  is  possible  or  not  to  get  all  the  eigenvalues  equal  so  as 
to  maximize  the  determinant  of  the  FIM.  If  not,  it  requires 
that  one  sensor,  which  is  much  closer  to  target  than  the  rest, 
be  n/2  to  all  other  sensors  that  must  stay  collinear. 
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Brining  (31b)  to  (29a)  gives: 


The  condition  for  the  remaining  sensors  to  be  collinear  is 
(31a)  equivalent  to  a  single  sensor  having  the  same  LOS  vector 
with  an  equivalent  range  to  target.  It  is  also  equivalent  to  a 
(31b)  single  sensor  making  multiple  independent  observations  so 
that  the  resulting  variance-weighted  inverse  range  is 
(31c)  comparable  to  the  closer  sensor.  Either  equivalence  enables 
FIM  to  be  diagonalized.  The  condition  becomes  a  two¬ 
sensor  scenario,  which  is  similar  to  what  is  presented  in 
Section  3.1  and  3.3. 
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A  generalization  of  the  condition  (43)  of  Theorem  9  in  [1] 
is  as  follows.  For  m  sensors  with  variance  gft  i  =  1,  ...,  m. 

The  conditions  of  Theorem  1  can  be  achieved  for  k 
measurement  if 
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and  for  sensor  j  that  produces  the  smallest  measurement 
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The  above  analysis  shows  a  closed  form  solution  for  2-D 
optimal  placement  of  2  sensors.  However,  there  may  be  no 
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Note  the  first  condition  (35a)  puts  a  lower  bound  on  a 
possible  k  and  the  second  term  (35b)  put  the  upper  bound. 
For  two  sensors,  this  is  equivalent  to  (34). 


V0.7375 

x/0.7375 

0.6557  0 

x/0.2625 
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0  0.6557 

4.  Simulation  Examples 

In  this  section,  three  examples  are  used  to  illustrate  the 
optimal  placement  conditions  derived  in  previous  sections. 

Example  1.  When  the  prior  covariance  is  not  circular, 
placing  sensors  evenly  all  around  the  error  ellipse  is  not  an 
optimal  approach.  Similarly,  optimizing  the  angular 
orientation  of  a  sensor,  one  measurement  at  a  time,  is  not 
an  optimal  strategy,  either.  These  two  points  are  illustrated 
by  the  following  numerical  example  in  contrast  to  the 
optimal  condition  given  in  Section  3.2. 

Consider  the  prior  covariance  P0  =  [20  0;  0  1]  and  two 
sensor  updates  with  equal  effective  measurement  error 
<72  =  1 .  Following  the  optimal  single  sensor  update  strategy 
of  selecting  the  angular  orientation  that  aligns  with  the 
maximum  eigenvector,  the  updated  error  covariance  is 
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"0.9524  0  " 
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This  choice  of  angular  positions  for  the  two  sensors 
surrounds  the  target  and  the  resulting  performance 
measures  for  P  are  tp\  =  trace( P)  =  1.4524  and  dp\  =  det( P) 
=  0.4762.  The  corresponding  performance  measures  for  M 
are  tm\  =  trace{ M)  =  3.05  and  dni\  =  det( M)  =  2.30. 


From  (36),  it  is  easy  to  see  that  the  placement  angles  are  ft 
~  -02  =  30.8203°.  For  this  placement,  the  performance 
measures  for  P  are  tp->  =  trace{ P)  =  1.3115  and  dp2  =  del(P) 
=  0.4300.  The  corresponding  performance  measures  for  M 
are  tm2  =  trace( M)  =  3.05  and  dm2  =  det( M)  =  2.3256. 
Clearly,  under  the  same  measurement  quality  (tm  \  =  tm2  = 
3.05),  the  second  choice  is  better  (dm2  =  2.3256  >  dni\  = 
2.3). 

Fig.  1  shows  the  color-coded  surface  of  the  determinant  of 
the  updated  FIM  (14a)  as  a  function  of  for  9\  and  ft  for  this 
example  with  M0  =  [0.05  0;  0  1]  and  cf2  =1.  Dark  red 
areas  indicate  large  values  whereas  dark  blue  areas  indicate 
small  values. 

As  shown,  there  are  four  pairs  of  maximum  values.  The 
first  pair  occurs  when  ft  =  30.8203°  and  ft  =  149.1797° 
and  329.1797°,  respectively.  Note  that  the  two  values  of  ft 
are  off  by  180°,  thus  being  along  the  same  direction.  For 
the  second  pair,  the  peaks  appear  at  ft  =  149.1797°  and  ft 
=  30.8203°  and  210.8203°,  respectively.  Again  the  two 
values  of  ft  are  off  by  180°.  It  is  clear  that  the  value  of  ft 
for  the  second  pair  is  in  fact  that  of  ft  in  the  first  pair,  ft 
and  ft  switch  their  positions  in  the  two  pairs.  Finally,  the 
third  and  fourth  pairs  are  repeated  patterns  of  the  first  and 
second  pairs  with  180°  in  ft. 


Flowever,  it  is  possible  to  position  the  sensors  to  provide  a 
tighter  covariance  error.  Namely, 

MQ  =  diag([0.05,  1],a*  =a*  =  3,  k=  1,  max(det(M))  =  0.71667  MQ  =  diag([0.05,  1],a*  =  a*  =  1,  k=  1,  max(det(M))  =  £3256 


Fig.  1.  Determinant  of  Updated  FIM  vs.  6\  and  ^ 
for  Given  Mo  and  (Example  1) 

MQ  =  diag([0.05,  1],o^  =  2,0*  =  6,  k=  1,  max(det(M))  =  0.71667 


ftj,  deg 

Fig.  2.  Determinant  of  Updated  FIM  vs.  ft  and  ft 
for  Given  M0  and  g2  (Example  2,  k  =  1) 

MQ  =  diag([0.05,  1],  o*  =  2,o*  =  6,  k  =  £  max(det(M))  =  1.42 


Fig.  4.  Detenninant  of  Updated  FIM  vs.  ft  and  ft  Fig.  5.  Detenninant  of  Updated  FIM  vs.  ft  and  ft 
for  Given  M0  and  g2  (Example  3,  k  =  1)  for  Given  M0  and  g2  (Example  3,  k  =  2) 


MQ  =  diag([0.05,  1],o*  =o*  =  3,  k=  2,  max(det(M))  =  1.42 


Fig.  3.  Detenninant  of  Updated  FIM  vs.  0\  and  ^ 
for  Given  Mo  and  tf2  (Example  2,k  =  2) 

M0  =  diag([0.05,  1],o*  =  2 ,o*=  6,  k=  3,  max(det(M))  =  £325 
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Fig.  6.  Detenninant  of  Updated  FIM  vs.  ft  and  ft 
for  Given  M0  and  g2  (Example  3,  A:  =  3 
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Example  2.  In  the  second  example,  we  still  have  M0  = 
[0.05  0;  0  1]  but  a2  =3-  It  is  easy  to  verify  from  (18)  that 
there  is  no  single-look  optimal  placement.  The  color-coded 
surface  of  the  determinant  of  the  updated  FIM  as  a  function 
of  for  ft  and  ft  for  this  example  with  k  =  1  is  shown  in  Fig. 
2.  The  largest  values  occur  at  9\  =  0°  (ft  =  180°)  and  92  = 
0°  and  180°,  respectively,  meaning  along  the  largest 
eigenvalue  direction. 

However,  if  the  sensors  take  two  independent 
measurements  per  placement,  the  resulting  updated  FIM  is 
shown  in  Fig.  3.  The  first  pair  of  peaks  occur  at  (ft  = 
22.2807°,  ft  =  337.7193°)  and  (ft  =  22.2807°,  ft  = 
157.7193°)  where  the  two  values  ft  are  off  by  180°, 
pointing  to  the  same  direction. 

The  second  pair  of  peaks  occur  when  ft  and  ft  just  switch 
their  values,  that  is,  (ft  =  157.7193°,  ft  =  202.2807°)  and 
(ft  =  157.7193°,  ft  =  22.2807°).  Again,  the  two  values  ft 
are  off  by  180°,  pointing  to  the  same  direction. 

The  third  and  fourth  pairs  repeat  the  first  and  second  pairs 
every  180°  as  discussed  in  Example  1  and  they  actually 
point  to  the  same  direction. 

Example  3.  In  the  third  example,  we  still  have  M0  =  [0.05 
0;  0  1]  but  a*2  =  2  and  &l=6-  It  is  easy  to  verify  from  (34) 
that  there  is  no  single-look  optimal  placement.  The  color- 
coded  surface  of  the  determinant  of  the  updated  FIM  as  a 
function  of  for  ft  and  ft  for  this  example  with  k  =  1  is 
shown  in  Fig.  4.  As  in  Example  2  (Fig.  2),  when  the 
optimal  condition  is  not  met,  the  peaks  occur  at  ft  =  0°  ( ft 
=  180°)  and  ft  =  0°  and  180°,  respectively,  meaning  along 
the  largest  eigenvalue  direction. 

However,  if  the  sensors  take  two  independent 
measurements  per  placement,  the  resulting  updated  FIM  is 
shown  in  Fig.  5  for  k  =  2.  The  first  pair  of  peaks  occur  at 
(ft  =  43.7034°,  ft  =  -10.4142°)  and  (ft  =  43.7034°,  ft  = 
169.5858°)  where  again  the  two  values  ft  are  off  by  180°. 
The  second  to  fourth  pairs  of  peaks  appear  following  the 
same  patterns  as  discussed  in  previous  examples. 

When  k  >  3,  the  optimal  condition  (34)  does  not  hold. 
Indeed,  the  resulting  updated  FIM  with  three  independent 
measurements  is  shown  in  Fig.  6.  As  shown,  when  the 
optimal  condition  is  not  met,  the  peaks  occur  at  ft  =  90°  ( ft 
=  210°)  and  ft  =  0°  and  180°,  respectively.  This  indicates 
the  desired  updates  are  along  the  two  eigenvalue  directions. 

5.  Conclusions 

In  this  paper  we  derived  conditions  for  optimal  placement 
of  two  heterogeneous  sensors  in  tracking  of  uncertain 
targets  over  several  time  steps.  The  conditions  were  derived 
based  on  the  maximization  of  the  updated  Fisher 
information  matrix  from  an  arbitrary  prior  characterizing 
the  uncertainty  about  the  initial  location  of  a  target.  The 


two  heterogeneous  sensors  can  be  of  the  same  or  different 
types  such  as  ranging  sensors,  bearing-only  sensors,  or 
both.  The  sensors  can  also  offer  measurements  of  different 
qualities.  The  results  presented  in  this  paper  complement 
other  recently  published  work  on  instantaneous  optimal 
placement  that  considered  same  sensor  types  and  same 
measurement  qualities  for  targets  at  perfectly  known 
locations  or  with  the  target  location  uncertainty  averaged 
out  via  the  expected  value  of  the  determinant  of  the  Fisher 
information  matrix. 

Future  work  will  consider  the  cases  where  two  cooperative 
mobile  sensors  take  their  measurements  neither  at  the  same 
rate  nor  at  the  same  time  where  two  numbers  of 
independent  update  k\  and  k2  can  used  to  account  for  this 
asynchronous  nature  of  operations  with  distributed  sensors. 
Another  effort  is  to  extend  the  results  of  this  paper  to  m  >  2 
sensors.  In  the  spirit  of  [2],  one  may  take  the  expectation  of 
a  matrix  measure  of  choice  and  one  choice  is  the  trace  of 
the  inverse  of  M.  It  is  of  great  interest  to  bring  in  such 
aspects  as  placement  costs  and  risks  into  consideration, 
which  will  be  another  direction  of  our  future  study. 
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